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Porównanie mierników współliniowości 


A Comparison of the Measurements of Colinearity 


Aktualny stan wiedzy na temat współliniowości zmiennych objaśniają- 
cych jest w literaturze bardzo obszerny. Szczególnie dużo jest publikacji 
na temat ujemnych skutków współliniowości w procesie estymacji param- 
terów strukturalnych metodą najmniejszych kwadratów. Do najważniej- 
szych prac dotyczących tego problemu można zaliczyć [4], [5], [11], [19]. 
Ponieważ zjawisko współliniowości ma szczególnie negatywny wpływ na 
szacowanie modelu klasycznymi metodami estymacji, dlatego w literatu- 
rze pojawiło się wiele prac na temat metod eliminujących wspólliniowość 
zmiennych. Sa to między innymi metoda regresji grzbietowej podana przez 
E. Hoerla i W. Kennarda w pracy [9], metody ortogoalizacyjne (por. prace: 
[7], [10], [12], [13], [16], [17], [18]), metoda różniczki zupełnej, opisana np. 
w pracy [14] i metoda transformacji zmiennych (por. [16]). Ze względu na 
brak precyzyjności w definicji współliniowości zmiennych objaśniających, 
która mówi, że zmienne objaśniające są wspólliniowe, jeżeli co najmniej 
jedna z nich jest silnie skorelowana z pozostalymi, opracowano w litera- 
turze dużą liczbę kryteriów, na podstawie których mierzy się natężenie 
współliniowości. 

Niniejszy artykuł poświęcony jest eksperymentowi symulacyjnemu do- 
tyczącemu porównań charakterystyk zmienności funkcji opisujących pewne 
wybrane mierniki wspólliniowości w zależności od natężenia współliniowo- 
ści między zmiennymi objaśniającymi liniowego modelu ekonometrycznego. 
Miernikami współliniowości, które badano w przeprowadzonym ekspery- 
mentcie byly:! 


1 We wszystkich rozważanych miernikach współliniowości k oznacza liczbę zmiennych 
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— integralny współczynnik zasobu informacji brutto 


k 


Q zz ys TĘX,» 


j=l 
gdzie ryx, jest współczynnikiem korelacji między zmienną objaśnioną 
Y a j-tą zmienną objaśniającą X;, 

— integralny współczynnik zasobu informacji netto 


k 
G= 90 gi 
j=l 


gdzie 
Bzryx, „gdy Bjryx, > 0 

0 „gdy Śzryx, <0’ 

przy czym B; jest oceną j-tego parametru strukturalnego 8; modelu, 
— integralny współczynnik zasobu informacji powtórzonej 


eat 
E 


F'=Q-G, 
— integralny współczynnik pojemności 
k 
Bt} hi, 
j=l 
gdzie 
yx, 


ne M 
E , 
ior pare eos Irx,x;| 


— statystyka Farrara—Glaubera 


2 
AP =" 


s 1, 
n-1- glk + 5) log det R, 
gdzie R jest macierzą współczynników korelacji między zmiennymi ob- 
jaśniającymi, zaś n oznacza liczbę obserwacji zmiennych modelu, 
— statystyka Haitovskiego 
1 
X = — [n -1- 5 (2k + 5) log(1 — det R), 


objaśniających modelu. Bliższe informacje na temat tych mierników można znaleźć np. 
w pracach: [1], [3], [5], [6], [8], [15]. 


Porównanie mierników współliniowości 
integralny współczynnik współliniowości 
z 2 
A=Q-fR'*, 


gdzie R? jest współczynnikiem determinacji, 
przeciętny współczynnik indywidualnej wspólliniowości 


Ai 
6 = (= Ji ) i 
= k ; 


gdzie 


det R 
êj z p z det R; 
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przy czym det Rj; jest minorem głównym macierzy R otrzymanym przez 


skreślenie j-tego wiersza oraz j-tej kolumny, 
przeciętny wspólczynnik Farrara—Glaubera 


gdzie 


przy czym r” jest j-tym elementem glównej przekątnej macierzy R`}, 


przeciętny wspólczynnik zasobu informacji brutto 


Eijatiti Tux, 
k(k = 1) 


? 
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— przeciętna wartość statystyk Studenta 


46) = = (455)? 


. 
k 


— przeciętny współczynnik regresji standaryzowanej 


PR z WZICJĄ 
"AW Td" 


gdzie B; jest oceną j-tego parametru strukturalnego 8; modelu o zmien- 
nych standaryzowanych na 0 — 1, 


— przeciętny wspólczynnik określoności indywidualnej 


- |= Dis (dj)? 
i- JEF 


gdzie 
d, A Â; Lisl: — Y)(Tij — Zj) 
? 
pai (Yi = 9)? 
przy czym y; , t;; są i-tymi realizacjami odpowiednio zmiennych Y oraz 
Xj, zaś 


o 1 n R 1 n 
ż;= n te J= „ŻW 


— przeciętny wspólczynnik zasobu informacji netto 


gdzie 
m; = = Bibjrx;x, , 


s=1 
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— przeciętny współczynnik integralnej koincydencji 


gdzie 


k 
Pi = 2 rx,x,did; = 2d, TY A; 
i=1 
zaś d; jest i-ta współrzędną wektora d = U”! Ro, przy czym Ro oznacza 
k-elementowy wektor o współrzędnych ry x,, natomiast U jest macierzą 
o elementach równych 


wj=ryx,ryx, (j= 1,2,:**,k;i£ j) 
uj =1 (j= 526,75 


— przeciętny wspólczynnik pojemności 


"WIR 
ka 
Aby móc przeprowadzić eksperyment symulacyjny na wyżej przedsta- 

wiony temat należało przede wszystkim zaprojektować taki model ekono- 

metryczny, w którym można by było zmieniać natężenie współliniowości 

w zbiorze zmiennych objaśniających. Modelem takim był model, dla któ- 

rego macierz korelacji między zmiennymi objaśniającymi miała następującą 

postać:? 


£ 
1 0 a *:: a 
0 1 0 a a 
R(k.p)(2) = 0 0 see I a ots @ , (1) 
a a *. a 1 © a 
i a a =œ- a a *:* 1 


„| kxk 
gdzie a = ry_pyij dla j £ k -—p+i(t=1,2,---,p; J = 1,2,:::,k). 
Natężenie współliniowości w zbiorze zmiennych objaśniających było 
zmieniane na skutek zmian wartości współczynnika a, tak by dla odpo- 
wiedniego a macierz R(, „)(a) była macierzą korelacji, czyli każdy jej minor 
główny musi mieć wartość większą lub równą zero. 


2 Taka lokalizacja p-elementowego zbioru wierszy została wybrana dla jasności prowa- 
dzonego wywodu, nie zmienia to jednak ogólności rozważań. 
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Aby wykazać jaki jest największy przedział określoności współczynnika 
a, by macierz (1) były macierzą korelacji sformulowano i udowodniono 
następujące dwa twierdzenia, które są potrzebne do tego celu. 
Twierdzenie 1 
Jeżeli k,p € N (N — zbiór liczb naturalnych) oraz 2 < p< k — 1, to 


det Ry,py(a) = (1-4)?! |-p(k- pja? +(p- Da+1], | (2) 


gdzie R(,,)(a) jest macierzą określoną wzorem (1). 
Dowód 

Aby wykazać równość (2) postępujemy według następującego algorytmu: 
(1) W wyznaczniku macierzy (1) od kolumny o numerze k—p+i odejmujemy 
kolumnę o numerze k — p + i + 1 dla każdego wskaźnika ż = 1,2,:::,p— 1. 
W wyniku tego przekształcenia otrzymujemy wyznacznik postaci: 


1 30 eee 0 0 0 Ke IG 
0 1 .:- 0 0 0 ©. a 
0 0 s 1 0 0 .. a 

+ a l-a 0 . a 
a a :':': a a—l l=a ::: a 


a a a wa 0 0 a 1 
(2) W powyższym wyznaczniku k — p pierwszych kolumn mnożymy przez 


—a i dodajemy kolejno do ostatniej kolumny. Po tym przeksztalceniu otrzy- 
mujemy następujący wyznacznik 


bo 0. 270 0 DL eee 0 
Ko <a M 0 WZ 0 
iene BI 0 GP 4. 0 
a a = a l~a 0  ::: a-(k-p)a? 
a a >- a a-l l-a «++ a—(k— p)a? 
Gat 337 awm 0 0 +++ 1-(k- pa? 


(3) W wyznaczniku otrzymanym powyżej kolumnę o numerze k- p + i 
mnożymy przez wyrażenie 


i(a — (k — p)a?) 
a 


p — (az 1) 
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dla i = 1,2,---, p— 1 i odejmujemy kolejno od k-tej kolumny. (Dla a = 1 
macierz (1) ma p jednakowych kolumn i wtedy jej wyznacznik jest równy 0). 
W wyniku tego przekształcenia otrzymujemy wyznacznik postaci 


i 0a: 0 0 0 z 0 
0 1 = 0 0 0 see 0 
0 0 1 0 0 ... 0 
a a a l-a 0 ee 0 
a a »- a a—l l-a œ- 0 


a a *':: a 0 0 +++ 1-(k-p)a + (p— 1)fa — (k — p)a?] 


Ponieważ jest to wyznacznik macierzy trójkątnej dolnej, więc jest on równy 
następującemu wyrażeniu 


(1 - a)P™ [1 — (k — p)a? + (p— 1)(a — (k — p)a?)]. 
Pokazaliśmy więc, że 
det Rixp(4) = (1 — a)? — (k — p)a? + (p — 1)(a — (k — p)a?”)] 
czyli ostatecznie 
det Fq,y(a) = (1 — a?7'|--pfk — p)a? + (p — 1)a + 1] 


c.b.d.o. 
Zanim przedstawimy drugie z zapowiedzianych twierdzeń udowodnimy 
następujący lemat, z którego będziemy korzystać w trakcie dowodu tego 
twierdzenia. 
Lemat 1 
Jeżeli pke N,2< p< k-l1 oraza € (-1,1), to: 


(1) ao € (—1,0)i bo € (0,1) (3) 
(2) det Ręp(aj=zO0$a=agVa=lVa=b (4) 
(3) det Ri, p)(a) > 0 & a € (ao, bo) (5) 
(4) det Ryta) <OSaeE (-1, ag) V (bo, 1), (6) 
gdzie 
NE Dała 1)? + 4p(k — p) fet Pls a V(p — 1)? + 4p(k — p) 
j 2p(k — p) 2p(k — p) 


(7) 
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Dowód 

Warunek (3) jest oczywisty na mocy elementarnych przekształceń. Nato- 
miast warunek (4) wynika wprost ze wzorów (7) oraz (2). Z kolei z zależności 
(2) oraz (4) otrzymujemy następujące równoważności: 


| a<lAag<a<bo 
det Rę.„y(a) > 0 © + a>lAao<a<boAp-l jest parzysta 
| a>lA(a<agVa>bo)Ap-1l_ jest nieparzysta 


oraz 


| a< 1A(a< agVa > bo) 
det Ry. p)(@) <0 4 a>1lA(a<aoAa>bo)Ap-1 jest parzysta 
| a>lA(ag<a<by)Ap-1 jest nieparzysta 


Stąd po uwzględnieniu warunku (3) otrzymujemy, że 

det Rę.p)( 4) > 0 & ao < c < bo V (a > LAp— L jest liczbą nieparzystą) 
oraz 
det Rix p (a) <0%a<aoVbo<a<1(a>1Ap- l jest liczbą parzystą) 


czyli równoważności (5) i (6) dla a € (—1,1) są prawdziwe, 

c.b.d.o. 
Twierdzenie 2 
Jeżeli pk e NI2< p< k—1, to przedział (ag,bo) jest największym 
przedziałem zawartym w (—1, 1) takim, że 


N det Rę_;7-3(a) > 0 (8) 
a€(ao,bo) 1=0,1,---,p—1 
gdzie ag, by określone są wzorem (7), natomiast Rq,_;„_:)(a) jest podmacie- 
rzą macierzy (1), która powstaje przez skreślenie w niej i-ostatnich wierszy 
oraz i-ostatnich kolumn. 
Dowód 
Niech a € (—1,1), wówczas z lematu 1 wynika, że 


N GE (aj, bi) > det Rę_;9-3(a) > 0, (9) 
t=0,--,p—1 
gdzie 
_p-i-l-y(p-i-1)*+4(p- i)(k — p) 
2(p — tj(k — p) 


a; 


(1=0,1,:::,p= 1) (10) 
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po 14 Vip= SEFA- ES p) (i 
2(p — i)(k — p) 

Z warunku (9) wynika, że aby zakończyć dowód wystarczy wykazać, że ciąg 

(a;) jest malejący, zaś ciąg (b;) niemalejący, gdyż wtedy 


b; = =0,1,::-,p= 1) (11) 


czyli 


N (ai, 0) = (ao, 0) oraz () (0,b;) = (0, bo) 


i=0,1,--,p—1 i=0,1,---,p-1 


Poniżej najpierw wykażemy, że ciąg (a;i) jest malejący czyli, że 


A  «w> Gina (12) 
i=0,1,++,p—2 


Ponieważ dla a € (—1,0) w myśl lematu I mamy warunek 
det Ri,-i-1,p-i-1)(@) > 0 $ a E (041,0) (i = 0,1,---, p — 2) (13) 
więc aby udowodnić nierówność (12) wystarczy wykazać, że 
det R(k_i-1,p-i-1)(56:) > 0 (1 = 0,1,--+,p — 2) (14) 


gdyż  € (—1,0) (na mocy lematu 1), czyli wtedy z równoważności (13) 
otrzymamy, że a; > aj4; (1=0,1,:::,p— 2). 

W celu wykazania nierówności (14) przeksztalcamy det Ry, 
stępujący sposób: 


~ip-i)(@) W na- 
det Rg 05 (6) u 
=(1-a)?*H-(p-i-1+1)(k — p)a? +(p-1—2+1)a+1] = 
=(1- a)P*"H-(p-i- 1)(k — p)a? + (p— i- 2)a+ 1] + 
+(1 - a?7""H-(k — pja? + a] = 
= (1 — a) det Rqe_j-ą,p-i-1(2) + (1 - a)77*7*[-(k — p)a? + a]. 
Stąd dla a ź 1 


det Rik-ip-i)(0) 


Toe 0a paka] (15) 


det Riepe) = 
Ponieważ na mocy lematu 1, a; ź 1 oraz det Rix-i,p-i)(ai) = 0, więc 


det R(p_;—1,p-i-1)(0:) = —(1- aj)? (k — pja? + aj]. (16) 
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Z faktu, że a; € (—1,0) oraz równości (16) wynika, że 
det Rix—i-1,p-i-1)(@i) > 0 & 1 > (k — p)a;, 
a stąd korzystając ze wzoru (10) otrzymujemy równoważność 
det Rqk_;-1,p-t-1)(8) ŻO%p<k-1(1=0,1,---,p— 2), 


czyli nierówność (14) jest prawdziwa. 
Z kolei udowodnimy, że ciąg (b;) jest niemalejący czyli, że 


A bi S dint (17) 


2=0,1,---,p~—2 

Poniewaz jak latwo sprawdzi¢ zachodzi nastepujacy warunek 

b= 16 p=k-—1(?=0,1,---,p=2) (18) 
więc gdy p = k — 1, to ciąg (b;) jest ciągiem stałym. 
Z lematu 1 wynika, że b, € (0,1) więc należy jeszcze rozważyć przypadek 
gdy b; € (0,1). Poniżej pokażemy, że w tym przypadku ciąg (b;) jest ciągiem 
rosnącym. 
Ponieważ dla b; € (0, 1) z warunku (5) wynika, że 

det Rę=i=1p1-1)(0:) > 0 eb; € (0,b;41) (1 = 0,1,---,p— 2), 
więc aby udowodnić nierówność 
bi < bign (1=0,1,:::,p=2) 
wystarczy wykazać następującą równoważność 
det Rik-i-1,p-i-1)(b:) > 0 & p<k-1(1=0,1,::-,p— 2). 


Zauważmy, że dla b; # 1 na mocy zależności (4) i (15) zachodzi równość 


det Rik-iz1p-i-1)(bi) = —(1 = biPTH- (k = p)b? + bi] (i = 0,1,-+-, p 2). 
(19) 
Z kolei wykorzystując wzór (11) i (19) otrzymujemy po przeksztalceniu, ze 


det Rgt Eb 0 Gp<k-L, 


c.b.d.o. 
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Na podstawie twierdzenia 2 widzimy, że macierz (1) jest macierzą 
korelacji wtedy i tylko wtedy gdy a € (ag,bo) lub a = 1. Natężenie 
wspólliniowości w zbiorze zmiennych objaśniających o macierzy korelacji 
określonej wzorem (1) możemy zmieniać w sposób regularny np. dzieląc ag 
lub bę przez dowolną liczbę naturalną m. Wtedy otrzymujemy następujące 
wzory: 


() _P-1-v(p- IP +4PlE=P),._y pao... mah) ( 

ait} — Pipl =o) (r—1) (r=1,2,---,m+1) (20) 
=n — ‘ p — 

pel B= PSEA OE (7 SETI WA CES (21) 


2mp(k — p) 
określające m + 1 kolejnych wartości współczynnika a odpowiednio dla 
a € (ag,0) oraz (bo, 0). 

Aby się przekonać jaka jest monotoniczność funkcji det Ri, »)(a) w prze- 
działe < ag,bo > udowodnimy następujące 
Twierdzenie 3 
Jeżeli p, k € N oraz 2 < p < k-1, to funkcja det R(x „)(a) rośnie w przedziale 
(ag, 0) oraz maleje w przedziale (0, bo). 
Dowód 
Korzystamy z następujących implikacji 


(1) Nzelznza) f'(x) < 0 = funkcja f jest malejąca w (21,22) (22) 
(2) Azete, z) /'(7) > 0 > funkcja f jest malejąca w (21,22) (23) 


Zauważmy, że pochodna funkcji det Rq.„y(a) liczona względem zmiennej a 
wyraża się wzorem 


(det Rę.(a)) =(1- a)?7?pl(k — p)(p + 1)a? — (2(k — p) + p — 1)a] (24) 
Z równości (24) wynika, że 


a<1A0<a< GETRY 


(det Rix p (a) <0 $ | a>1AQ0<a< LS A p jest parzystaV 


a>lA(a<0Va> AZ 


G-p(p+1) A p jest nieparzysta 


Ale dla p < k — 1 prawdziwa jest nierówność 


p-—1+2(k — p) 


"<b per =" 


(25) 
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wiec 


(k=p)(p+1) 


| a> lAp jest liczbą nieparzystą 


| O<a< ZER 
(det Rq.y(a)) <0 © 


Ponieważ jak łatwo sprawdzić zachodzi nierówność 


p-l+v(p-1)*+dplk-p) „p-1+2(k— p) 
2p(k — p) “(k= p)(p+1) ° 


więc z warunku (22), (26) oraz powyższej nierówności wynika, że jeśli 


= aq v2 = 
eml 1+ J(p— 1)? + 4p(k DIY 
2p(k — p) 


to det Ri, p)(a) jest funkcją malejącą. 
Z kolei 


a<0 
V 
(det Rę.y(a)) > 0 & 4 p-1+2(k-p) 
| <a <a 
| a> lAp jest liczbą parzystą, 


czyli z warunku (23) wynika, że jeśli 


po b= Ya + (FE =D) 0) 
2p(k — p) "Z 


to det Ry;,„)(a) jest funkcją rosnącą. 


a E< 


c.b.c.o. 
Z twierdzenia 3 wynika, że zmieniając wartości współczynnika a w myśl 
wzorów (20) lub (21) dla kolejnych wartości r powodujemy, że wartość wy- 
znacznika macierzy (1) zmienia się w sposób malejący od maksymalnej 
wartości równej 1 do minimalnej równej 0, tzn. natężenie wspólliniowości 
zmienia się do najmniejszej do najwiekszej osiągalnej wartości. W celu wy- 
generowania k-wymiarowej zmiennej losowej (k zmiennych objaśnia jących) 
o rozkładzie normalnym z wartością oczekiwaną p = (41, /t2, * * * „jt ) i odchy- 
leniem standardowym o = (01, 02,°-+,0%) oraz macierzą korelacji Rix,p)(@) 
korzystano z następującego wzoru? 


X=CY +h, 


3 Por. pracę [2]. 
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gdzie Y = (Vy,Y2,::*,X4) jest k-wymiarową zmienną losową, dla której 
Y;() = 1,2,:-:,k) są niezależne i mają jednakowy rozklad N(0,1) zaś C 
jest macierzą nieosobliwą taką, że CCT jest macierzą wariancji-kowariancji 
zmiennej losowej X. 

Dla każdej wartości a()(r = 1,2,---,m + 1) wstawionej w miejsce a 
w macierzy (1) otrzymywano, wykorzystując powyższy generator k-zmien- 
nych objaśniających oraz zmienną objaśnianą a następnie odpowiadający 
im liniowy model (takie rozważano w niniejszej pracy) estymowano metodą 
najmniejszych kwadratów. Po estymacji obliczano wartość wszystkich ba- 
danych mierników wspolliniowosci. Ponieważ opisane postępowanie powta- 
rzano m+1 razy więc dla każdego miernika wspólliniowości otrzymano m+ 1 
różnych jego wartości, odpowiadających kolejno coraz większemu natężeniu 
wspólliniowości. Aby otrzymaną w ten sposób zależność danego miernika 
wspólliniowości od natężenia współliniowości opisać odpowiednią funkcją 
natężenie to wyrażano za pomocą wyznacznika macierzy korelacji. Następnie 
dla każdego miernika współliniowości aproksymowano otrzymane dla niego 
wartości wykresem tej funkcji wybranej spośród kilkunastu różnych typów 
funkcji, dla której współczynnik indeterminacji był najmniejszy. Potem dla 
wszystkich analizowanych mierników wspólliniowości obliczano wartości na- 
stępujących charakterystyk zmienności funkcji: 

— średniej prędkości względnej vsw oraz średniego przyspieszenia 
względnego asw zmiany wartości funkcji w przedziale (0,1), 

— prędkości względnej vw oraz przyspieszenia względnego ayy zmiany 
wartości funkcji w punkcie 0 i 1, 

— średniej elastyczności es funkcji w przedziale (0, 1) oraz elastyczności 
e funkcji w punkcie 0 i 1. 

Jeżeli chodzi o uporządkowanie mierników wspólliniowości wedlug miary 
lvsw|, to ukladają się one w następującej kolejności (od najmniejszej do 
największej) 


H,H,Q,6,t(8),G,6,8.Q,d, m, F', 9, Vy. P, F, À, Vw, 


przy czym wartości miary |vsw| dla mierników A, x% oraz © dość znacznie 
różnią się od wartości tej miary dla pozostałych mierników. Ze względu 
na miarę |vw(1)| rangi od najmniejszej do największej otrzymują kolejno 
mierniki: 


Q.6,:(8),d, m, 8,8, H, Ä,Q,G,F',g,Đ, X3 F, Xh, GA. 


Z kolei w wyniku porządkowania mierników wspólliniowości ze względu na 
bezwzględną wartość miary vw(0), rangi od najmniejszej do największej 
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otrzymuja kolejno mierniki: 
H, H,Q,6,G,5,Q, F',G,t(5),d, m, 5, r, .xF+9,0, Xp. . 


Miara |vw(0)| dla mierników x$, g, ©, Xy przyjmuje zdecydowanie większe 
wartości aniżeli w przypadku pozostałych mierników współliniowości. Jeżeli 
chodzi o miarę elastyczności średniej oraz elastyczności w punktach 1,0 
to mierniki wspólliniowości otrzymują analogicznie rangi jak odpowiednio 
w przypadku miar vs, tw(1), vzy(0). Z kolei miary |asw|, |aw(1)|, |aw(0)| 
przyjmują największą wartość dla miernika w, natomiast najmniejszą dla 
miernika H. Uporządkowanie wszystkich mierników współliniowości pod 
względem wyżej wymienionych miar przedstawia tablica 1. 


Tab. 1. Bezwzględne wartości miar asw, aw(1), aw(0) dla mierników współliniowości 


The absolute value of measurements asw, aw(1), aw(0) for the measurements of 
collinearity 


[owt | lew (OF 
0,196 0,573 
0,196 0,573 
0,573 0,573 
0,573 „0573 
1,640 1,099 
2,244 1,146 
2,244 1,146 
2,244 1,146 
2,244 1,146 
2,244 1,146 
2,244 1,146 
2,244 1,146 
2,982 1,283 
5,716 2,085 
5,716 2,085 
5,891 2,124 
9,197 3,037 

91,847 7,041 

802,158 21,774 


Ionn 


So a 


B 
Q 
F' 
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~ 
= 
= 


Ewa O 


Źródło: opracowania własne. 


Z przeprowadzonych badań wynika, że mierniki współliniowości reagują 
na zmianę natężenia wspólliniowości bardzo różnie. Największą „zmienność” 
wartości wskutek rosnącego natężenia wspólliniowości wykazują mierniki w, 
VF, natomiast najmniejszą (oprócz miary vyy(1)) mierniki H, H, 6. 


(1) 


[2] 
[3] 


[4] 
[5] 


[6] 
[7] 
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SUMMARY 


The present paper is devoted to a stimulating experiment concerning a comparison of 


the characteristics of variability of functions describing certain selected measurements of 
collinearity between the interpretative variables of the model. 


